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We propose a class of three-dimensional metamaterial architectures composed of a single doped semiconductor 
(e.g., n-Si) in air or vacuum that lead to unusual effective behavior of the classical Hall effect. Using an 
anisotropic structure, we numerically demonstrate a Hall voltage that is parallel—rather than orthogonal—to 
the external static magnetic-field vector (“parallel Hall effect”). The sign of this parallel Hall voltage can 
be determined by a structure parameter. Together with the previously demonstrated positive or negative 
orthogonal Hall voltage, we demonstrate four different sign combinations. 


In the textbook version of the classical Hall effect 1 , 
the Lorentz force due to an electric current in a semicon¬ 
ductor or metal plate and a magnetic-induction vector 
normal to the plate lead to an induced voltage. This 
Hall voltage is orthogonal to the current direction and 
orthogonal to the magnetic-field vector. In the linear 
regime, the Hall voltage is proportional to the magnetic 
field and proportional to the electric current. In isotropic 
bulk materials, the sign of the Hall voltage is determined 
by the dominant type of charged carrier, i.e., it is nega¬ 
tive for n-doped semiconductors (or free-electron metals) 
and positive for p-doped semiconductors. 

Recently, building upon the pioneering work of Bri- 
ane and Milton 2 , we have shown 3 that the sign of the 
Hall voltage can be modified by structure in a three- 
dimensional simple-cubic metamaterial lattice composed 
of intertwined tori made of only a single constituent ma¬ 
terial (e.g., n-doped silicon) in vacuum/air. We have 
argued that structures like that should be accessible 
experimentally by using current state-of-the-art three- 
dimensional fabrication technology. Due to the simple- 
cubic symmetry, the electrical properties have been 
isotropic 2,3 . 

In this Letter, we show that an anisotropic single¬ 
component structure can lead to a Hall voltage paral¬ 
lel (rather than orthogonal) to the magnetic-field vector. 
This finding complements the existing versions of the Hall 
effect, including the classical sign-inverted Hall effect 2,3 , 
the classical dynamic Hall effect 4 , the classical spin Hall 
effect 5 , the quantum Hall effect 6 , the fractional quantum 
Hall effect 7 , the quantum spin Hall effect 8 , and the opti¬ 
cal spin Hall effect 9 . 

The mathematical treatment underlying this Letter 
has been described in Ref. 3 and shall thus be sum¬ 
marized only briefly here: Using the software package 
COMSOL Multiphysics (MUMPS solver, more than 10 4 
tetrahedral elements per unit cell and a relative toler¬ 
ance of 10 -4 ), we numerically solve the static continu¬ 
ity equation V • = 0 for the electrostatic poten¬ 

tial <f> = (j>(r) with the conductivity tensor <r = a(r). 
In the linear regime and for a magnetic-induction vector 


B = (0, 0, B z \ it is given by 
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The current I x enters via the current density 
j = cfE = —crS/(f) through the boundary condition 
I x = — f A j • d A, where the integration plane A cuts the 
sample parallel to the yz- plane. To be concrete, we use 
material parameters corresponding to n-doped silicon as 
in 3 with scalar conductivity <to = 200 AV _1 m _1 and Hall 
coefficient Rh = —624 • 10 -6 m 3 A -1 s -1 , and a = 0 for 
the air voids. However, our results can easily be scaled 
to other parameter sets. Our calculations are performed 
for a metamaterial Hall bar composed of a finite num¬ 
ber of metamaterial unit cells N x , N y , and N z in the 
three spatial directions. Together with the cubic lattice 
constant a, this leads to the spatial extents of the Hall 
bar L x — aN x , L y = aN y , and L z = aN z . Large metal 
contacts (precisely, equipotential surfaces) parallel to the 
yz- plane covering the entire metamaterial faces are ap¬ 
plied to impose the current I x . For clarity, in this Let¬ 
ter, no metal contacts are explicitly added in the calcu¬ 
lations to pick up the Hall voltages. We have previously 
shown 3 that the corresponding perturbations are small. 
The shown Hall voltages are picked up at the unit cell 
in the center of the corresponding metamaterial sample 
surface. 

For the static case discussed in this Letter, the elec¬ 
tromagnetic wavelength, A, is formally infinitely large 
and, hence, an effective-medium description requiring the 
condition A/a >> 1 to be fulfilled is meaningful for any 
lattice constant a, e.g., for micrometer and centimeter 
scales. However, macroscopic structures are less favor¬ 
able because the Hall voltage is proportional to the Hall 
resistance, which, like any resistance, scales inversely pro¬ 
portional to size (if all dimensions are scaled simultane¬ 
ously). This means that, to obtain absolute Hall volt¬ 
ages that are readily measurable, macroscopic structures 
require orders of magnitude larger absolute currents I x 
and/or larger magnetic fields B z . Furthermore, the mag¬ 
netic fields obviously need to be homogeneous over much 
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FIG. 1. (a) Illustration of the suggested three-dimensional 
anisotropic metamaterial exhibiting a Hall voltage U z paral¬ 
lel to the magnetic-field vector (0,0, B z ) in addition to the 
usual orthogonal Hall voltage U y . The metamaterial is com¬ 
posed of a single constituent material, e.g., n-doped silicon, in 
air/vacuum. A constant current I x is imposed through the ar¬ 
ray of rods oriented along the x-direction. (b) One unit cell of 
the simple-cubic translational lattice with lattice constant a. 
The indicated parameters d y and d z can be negative/positive 
and determine the signs of U y and U z , respectively. The three 
different functional parts are made of a single material, they 
are colored for illustration only. The current-carrying straight 
rods are highlighted in red. The structures highlighted in 
blue and green lead to Hall voltages in the z- and ^-direction, 
respectively. The unit cell shown corresponds to a positive 
value of d z and a negative value of d y , the latter leading to 
an inversion of the orthogonal Hall voltage. 
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FIG. 2. Calculated Hall voltages for the metamaterial struc¬ 
ture illustrated in Fig. 1. (a) U z versus d z for three differ¬ 
ent values of d y . (b) U y versus d y for three different val¬ 
ues of d z . Parameters are: do — 200AV -1 m -1 , Ru = 
-624 • 10 -6 m 3 A -1 s -1 , I x = 0.1mA, B z = 1T, N x = 9, 
N y — N z = 3, h = 15 pm, r — 2.5 pm and a = 40 pm. 


larger regions in space for macroscopic structures. We 
have thus chosen a lattice constant of a = 40 pm in all of 
the calculations shown in this Letter. We believe that the 
corresponding three-dimensional structures can be fabri¬ 
cated by using three-dimensional direct laser writing (as 
discussed in more detail in Ref. 3). 

Intuitively speaking, the sign of a Hall voltage can sim¬ 
ply be flipped by interchanging the two wires used to pick 
up this Hall voltage from the Hall bar. To avoid short cir¬ 
cuits due to crossings of wires, this interchanging requires 
to go out of a common plane, i.e., it requires a non-planar 
three-dimensional geometry 2 . In the metamaterials pre¬ 
viously suggested in Ref. 3, this interchanging has rather 
been accomplished locally inside of the metamaterial by 
intertwined tori within each three-dimensional metama¬ 
terial unit cell. 

Expanding along this line, the novel blueprint depicted 
in Fig. 1 uses an array of straight semiconductor rods 
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d y = 8|jm dy = 0 |jm d y = -8 pm 



FIG. 3. Nine examples of Hall potential maps, i.e., electrostatic potential difference — 0(r, FL = 1 T) — 0(r, B z — 0) for 

N x = 9, N y = N z = 5 and different values of d y (different columns) and d z (different rows), plotted on a common false-color 
scale. The parameters are the same as in Fig. 2. The black arrows indicate the direction of the Hall-potential gradient. 


oriented along the x-direction (red in panel (b)). A to¬ 
tal electrical current I x flows via these rods through the 
entire sample. As a result, one locally gets an ordi¬ 
nary Hall voltage across each of these wires along the 
^/-direction. The structure contains two additional sets 
of “pick-up wires” made of the same material. One set of 


wires (green) connects the individual Hall voltages, which 
can be seen as local voltage sources, and adds them all 
up along the ^/-direction, leading to the usual orthog¬ 
onal Hall voltage U y . Likewise, the local voltages add 
up along the ^-direction through the second set of wires 
(blue), leading to the voltage U z parallel to the magnetic- 
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field direction. The two crucial geometry parameters are 
the quantities d y and d z (see Fig. 1). If they are nega¬ 
tive, the sign of each local voltage source is reversed and 
so are the macroscopic voltages U y and U z , respectively. 
The results of our numerical calculations shown in Fig. 2 
clearly confirm this conjecture. 

An illustrative way to see what is going on inside of 
the Hall bar is to inspect the Hall potential. We define 
the Hall potential as the difference of the electrostatic 
potential with magnetic field and that without. In this 
fashion, we subtract the overwhelmingly large but trivial 
potential ramp due to the applied bias voltage leading 
to the imposed current I x . The examples of the Hall- 
potential maps depicted in Fig. 3 demonstrate that the 
Hall voltages always build up from the “bulk” as op¬ 
posed to being a mere surface effect. Furthermore, the 
data in Fig. 3 show that we gain control of the orienta¬ 
tion of the Hall-potential gradient within the yz -plane 
(see black arrows). Equivalently, we can say that we 
have achieved unusual effective metamaterial conductiv¬ 
ity tensors a (see above) with adustable non-zero ele¬ 
ments <j xz = — G zx and a xy = —cr yx . 

As usual, the size and the aspect ratios of the Hall 
bar play an important role 10 regarding whether one ob¬ 
serves “bulk” effects or whether edge effects dominate. 
Figure 4 illustrates the scaling of the parallel Hall volt¬ 
age U z versus the length-to-height ratio of the Hall bar 
(i.e., versus N x /N z ). The depicted results resemble the 
typical behavior of the geometrical correction factor for 
ordinary bulk materials 10 . Additionally, it becomes clear 
that considering Hall bars composed of merely 9x3x3 
unit cells (as used for the calculations depicted in Fig. 2) 
is already sufficient to approach convergence towards case 
of infinitely many unit cells. 

For such sufficiently large and suitably shaped 
samples 10 and constant current I x , the ordinary orthog¬ 
onal Hall voltage (J7h = U y ) is proportional to l/L z and 
neither depends on L x nor on L y . In contrast, the paral¬ 
lel Hall voltage U z is proportional to 1 /L y and does not 
depend on L x and L z . Here, the increase of the Hall volt¬ 
age with the number of unit cells along the ^-direction 
is compensated by the fact that the current through one 
unit cell scales inversely proportional to the number of 
unit cells along z. The current through one unit cell and 
hence U z also scales inversely proportional to the number 
of unit cells along the ^/-direction. We have numerically 
verified all of these scaling laws (not depicted). 

In the special case of d y = 0 and d z = h (see Fig. 1), we 
obtain a metamaterial exhibiting the parallel Hall effect 
for an arbitrary orientation of the magnetic-field vector in 
the yz -plane (not depicted). In this special case, the class 
of structures we suggest shares some similarities to the 
structure suggested earlier by Briane and Milton 11 . How¬ 
ever, their structure requires two different constituent 
materials, one of which has anisotropic electrical prop¬ 
erties. In sharp contrast, our structure is composed of 
only a single locally electrically isotropic constituent ma¬ 
terial (like, e.g., n-doped silicon) - which is much simpler 
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FIG. 4. Scaling of the parallel Hall voltage versus the length- 
to-height ratio N x /N z . Parameters are N y = 1, d y = —8 pm 
and d z = 10 pm. The other parameters are as in Fig. 2. 

to realize experimentally. Furthermore, it is conceptu¬ 
ally even more striking to obtain the parallel Hall effect 
using only a porous structure of a single bulk constituent 
material. Finally, we note that the qualitative results de¬ 
scribed in this Letter can also be obtained by using hol¬ 
low rods instead of massive rods (not depicted). This fact 
enlarges the possibilities for experimental realizations. 

In conclusion, we have suggested three-dimensional 
metamaterials effectively exhibiting the parallel Hall ef¬ 
fect. These metamaterials are composed of only a single 
constituent material in vacuum/air, i.e., they can be con¬ 
sidered as a porous version of a bulk material. Broadly 
speaking, our results exemplify that the effective prop¬ 
erties of metamaterials can be way off those of the con¬ 
stituent (s). 
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